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Abstract. 

We develop a perturbative approach to study the supersymmetric non-linear sigma 
model characterized by a generic coupling matrix in the strong coupling limit. The 
method allows us to calculate explicitly the moments of the eigenfunctions and the 
two-level correlation function in the lowest order of the perturbative expansion. We 
find that the obtained expressions are equivalent to the results derived before for 
the corresponding random matrix ensembles. Such an equivalence is elucidated and 
generalized to all orders of the perturbative expansion by mapping the sigma model 
onto the field theory describing the almost diagonal random matrices. 
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1. Introduction 

One of the most powerful approaches to study complex quantum systems is based on 
the field-theoretical formulation. The non-linear cr-model, as an effective field theory 
describing long-distance physics, was first derived by Wegner [U |2] in the context of 
disordered quantum systems using the replica trick. An alternative formulation of the 
cr-model given by Efetov [3j 0] employed supersymmetry method making use of both 
commuting and anticommuting variables. The latter version of the cr-model can be 
successfully applied not only for problems, which can be solved on the perturbative 
level, but allows also for a non-perturbative treatment of the theory. Moreover, the 
supersymmetric cr-model can be derived as an effective field theory for various random 
matrix ensembles. In particular, it was used to prove that the statistical properties of the 
energy levels in small metallic disordered grains can be described by the Wigner-Dyson 
random matrix theory (RMT) [31 H] . 

The partition function of the generic supersymmetric non-linear cr-model on a lattice 

Z = J DQ e~T^"^tk U mhStiQ m Q k -J2 n So[Qn] ^ 

is defined in terms of supermatrices Q m associated with mth lattice site, a symmetric 
coupling matrix U m k > and a coupling constant t. S comprises all the local terms of 
the action including source terms. An explicit parametrization of Q and of the measure 



DQ = Yl dQ m is given in Appendix A 



An important feature of the cr-model is the fact that the super-matrix field Q 
satisfies a non-linear constraint Q 2 = 1. It is this constraint, which makes generally the 
analysis of the cr-model to be highly non-trivial. One possible way to tackle the problem 
is to use the perturbation theory of the non-linear cr-model in the metallic regime. The 
latter is characterized by a small value of the coupling constant t, which corresponds 
to high electrical conductivity in the case of disordered conductors. In this regime the 
statistical properties of energy levels and eigenfunctions are close to the corresponding 
ones of the Wigner-Dyson RMT. A large number of important and interesting results 
have been obtained following this route (see [5] and references therein). 

The opposite limit of the strong coupling has received much less attention. One 
of the reasons for that is the fact that in many cases the mapping between original 
microscopic models and the cr-model can be justified only at small values of the coupling 
constant. At the same time it is possible to construct microscopic models (such as 
for example n-orbital Wegner model [6]), which can be mapped onto the cr-model at 
all values of the coupling constant. This situation can be realized, for example, in a 
cluster of coupled disordered grains or chaotic systems, provided that each element of 
a cluster is well described by the Wigner-Dyson RMT. In such a case the properties 
of the cr-model in the weak and in the strong coupling limits are equally important. 
Moreover, some properties of the cr-model in the strong coupling limit might be of 
interest independently of whether or not it can be directly derived from a microscopic 
model. This is especially true for critical exponents of the cr-models which remain 
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critical at large values of the coupling constant, as critical exponents are the same for 
all microscopic models belonging to the same universality class [7]. 

Another reason, which makes the strong coupling limit to be less explored, is 
that technically this problem is more challenging, as the standard perturbation theory 
approach can not be applied there. This is in contrast to the weak coupling limit, where 
after separating out slow and fast modes, the latter can be integrated out order by order 
using the Wick's theorem. 

In this paper we develop a perturbation theory for the supersymmetric non-linear 
a- model defined on a lattice in the the strong coupling limit. The general idea of our 
approach is similar to the one of the virial expansion method introduced in Ref. [8] 
for studying almost- diagonal random matrices. After presenting a construction of 
the virial expansion for the generic cr-model, characterized by an arbitrary coupling 
matrix U, we calculate explicitly the first non-trivial terms of the expansion for two 
different quantities: the moments of the eigenf unctions I q and the two-level correlation 
function R^{u). Then we compare the obtained results with the corresponding 
expressions derived for almost-diagonal random matrices. Our main conclusion is that 
the expressions for I q and i?2^) are in one-to-one correspondence with the random 
matrix results. This fact can be traced back to the properties of the target space of Q 
and its counterpart discussed in Ref. [9]. We show that such a correspondence can be 
established at any order of the virial expansion under very general assumptions. 

The layout of this paper is as follows. In Sectional we present the general framework 
for the virial expansion formulated for the generic non-linear cr-model. Section [3] 
describes the calculations of the moments of the eigenfunctions I q and the distribution 
function of the eigenfunction intensities P(z) in the lowest non-trivial order of the 
coupling constant. Similar calculations for the two-level correlation function (^) 
are presented in Section 0J In Section [5] we apply the general formulas to the power- 
law non-linear cr-model and compare the obtained expression with the results known 
for the almost-diagonal random matrices. Finally, in Section [6] we discuss the origin of 
the results and demonstrate the equivalence between the cr-model description and the 
supersymmetric approach developed in Ref. [8] in the strong coupling limit. 

2. Virial expansion method 

The main idea of the virial expansion method is to represent the "interaction" part of 
the action j Ylim^k ^mk StrQ m Q fc as a sum over contributions characterized by a given 
number of "interacting" matrices Q m [8]. To this end, we denote by S[Q m ,Qk] = 
jUmk StrQmQk the sum of two terms which couple Q m with Q^. Then the "interaction" 
part of the partition function can be expanded as follows 



e -H m>k S[Q m ,Qu] = J| e ~S[Q m> Q k ] = JJ (1 + F mk ) = 



m>k m>k 




Virial expansion of the non-linear sigma model in the strong coupling limit 



4 



where F m k stands for e~ s ^ m '^ k ^ — 1. The first term of the above expansion is trivial. 
In the second one the contributions of pairs of supermatrices are collected together, the 
third one takes into account the contributions of triplets and so on. As a result we 
obtain the virial expansion of the partition function: 

Z= f DQe-ZnSo[Qn) + J2 j DQ ( e ~SKW^ - l) e -£ n So[Q„] + . . . .(3) 

m>k 

In the case of the disordered metallic grains the above expansion is useful, when the 
grains are weakly coupled. 

This general scheme is elaborated in the next two sections, where it is used to 
calculate the moments of the eigenfunctions I q and the two-level correlation function 
Ri{ui) in the unitary symmetry class ((3 = 2). 



3. Moments of the eigenfunctions 

The expression for the moments I q of the eigenfunctions (also known as generalized 
inverse participation ratios) in terms of the cr-model reads [TOl [TT] 

I q = lim^r f Str q (Q n kA) e-^^ UmpStrQmQp - r '^ StrAQm , (4) 

n 

where k = l<g)diag(l, —1) and A = diag(l, —1) in the retarded-advanced representation 
and N is the total number of sites. 

Applying the general scheme of the virial expansion outlined in the previous section, 
we obtain the following expansion for I q : 

Wf + + (5) 
The first term in this series corresponds to the first term of Eq. ([3]) and is given by 

J(°) = — J2 n m ^j- / D Q Str<? (Qn fcA) e~ v Em StrA «- . (6) 

One can notice that the integrand can be factorized and the integrals over different 
supermatrices Q m can be calculated separately. Moreover all the integrals over Q m with 
m 7^ n are equal to unity due to supersymmetry. The only non-trivial integral over Q n 
reproduces the standard expression for the moments of eigenvectors in the Wigner-Dyson 
RMT [11] . Thus we find that Iq ^ coincides with the standard RMT result 

Jf=r(g+1). (7) 

Physically it corresponds to the statistics of the eigenfunctions in a system of completely 
decoupled metallic grains. 

The second term of the expansion (jHJ) corresponding to the second term in Eq.([3]) 
reads 

/« = J2 Yl lh S IaT / DQ Str " ( QnkA ^ (e~ S[Qm ' Qp] ~ 1) e-" E '- StrAQm (8) 

n m>p 
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Again one can notice that, if m 7^ n and p 7^ n, then the double integral over Q m and 
Q p is determined solely by the anomalous contributions and hence equals to zero. Thus 
the terms with m = n or p = n are the only ones which should be taken into account: 

I« = V lim [ dQ m alQ n Str 9 (Q n kA) ( e - 5 K?-,Q«] - l) e -^trA(Q ra +Q n ) (g) 

^ — 4 r;— >0 4iV / 

m,n 

Now we need to calculate the double integral in the above expression. This 
task is straightforward, but technically quite challenging appreciating a non-trivial 
parametrization of the supermatrices. Fortunately, exactly the same integral appears 
at an intermediate step of the derivation of the non-linear a-model associated with the 
ensemble of the banded random matrices [TUl [12] ■ It turns out that the integrals over 
the Grassmann variables and the phases can be calculated exactly. The next step, which 
simplifies evaluation of the integral significantly, is to rescale the bosonic variables Xb 
in such a way that the limit rj — > can be taken explicitly in the integrand. Borrowing 
these results we arrive at the following expression for Iq 1 . 



IP 



" POO POO 

<7(<7 -1 )/ dz I duF mn (z,u) -T(q + 1) 

TO 'It t/ «/ 



1 

m,n 



sinli7 / 1\ sinh7 

cosh 7 h \u -\ — 

27 V u 



(10) 



[ID 



where parameter 7 = 4U mn /t is introduced in order to lighten the notation. The 
appearance of T(q + 1) in Eq. fllOl) is not accidental and can be traced back to the 
term containing —1 in Eq.()9]), which gives exactly the same contribution as the one 
calculated for Iq. 

The expansion (jHJ) makes sense only if its different terms scale as different powers of 
1/t. Therefore from the exact expression for 1^ derived above, we need to extract only 
the first non-trivial term in the limit t — > 00 or 7 — > 0. At this point our calculations 
depart from the a- model derivation, which assumes the opposite limit 7—7-00. 

It is obvious from the form of the function F mn that the limit 7 — > can not be 
taken before at least one of the integrations in Eq.( !T0|) is performed. One way to deal 
with this problem is to evaluate the integral over u exactly: 



duF mn (z,u) = e-v^V- 27 ^ + 7) CQSh7 + , ( * + lWl{2z + 7) Sinh7 .(12) 

7(22 + 7) 

The expansion of the last expression into the power-series in ^7 yields 

du F mn (z, u) = e- z z*- 2 [\ - ^ + 0(7)^ . (13) 
Substituting this result back into Eq. fllOl and integrating over z we obtain 

^ 1) = -^E\/|^- 1 ) r ^- 1 /2) + 0(7), q>\. (14) 
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The condition q > 1/2 is necessary for convergence of the integral at z — > 0. Restoring 
the dependence of 7 on U mn and taking into account the expression for Iq ^ we find the 
final result for the moments of the eigenfunctions 

\ mn / 

which is valid for an arbitrary coupling matrix U mn . 

The appearance of the factor T(q + 1) corresponding to the standard RMT result 
suggests that the eigenfunctions can be written as a product of two contributions. 
One of them describes the fluctuations of the eigenfunctions inside a grain and it is 
given by the Wigner-Dyson RMT. Another one, represented by the series in Eq.f fT5|) 
describes a change of the amplitudes of the wavefunctions at different grains. The same 
factorization occurs for the statistics of the eigenfunctions of one-dimensional and quasi- 
one- dimensional disordered wires and for arbitrary dimensionality in the metallic regime 

It follows from Eqs. (1711101113]) that the distribution function of the eigenfunction 
intensities is given by 



™=^( i -^4ev^+°<«- 1 ) )■ (16) 

\ mn / 

We would like to notice that the first two terms of the series presented in Eq. (|T5|) are 
in one-to-one correspondence with the expansion for the moments of the eigenfunctions 
found for almost diagonal random matrices characterized by an arbitrary variance matrix 
(\H mn \ 2 \) [IB] . It implies, in particular, the same expressions for the fractal dimensions 
in critical ensembles, as we discuss in Section [5] 

4. Two-level correlation function 

In this section we focus on the correlation function Rz(u)) measuring the correlations 
between two levels separated by an energy interval uj 

~ ^yp > ( 17 ) 

where v(E) is the density of states and the averaging over a random matrix ensemble 
or different disorder realizations is assumed. 

It is well known that ^(w) can be expressed as the following correlator of Q 
matrices [U [3] 

2 

' Str (Q n kA) 



-1 2 

#2M = ( 7^7 ] Re / DQ 



AN 



e -\ E m ^ p U mp StrQ m Q p +iip £,„ StrAQ m Qg-j 



where it is assumed that u has an infinitesimal imaginary part to ensure the convergence 
of the integral. This formula allows us to find the general form of the first two terms of 
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4V) 



i 

AN 



Re / DQ 



= (Jft) Re J DQ 



Str {Q n kA) 

str(g n M) 



e ^ Em strAQrl 



m>k 

In the following two subsections we evaluate (a;) and i2j 0^) respectively. 



(19) 
(20) 



4-1- Two-level correlation function in the zeroth order approximation 

Let us first calculate R^\uj). To this end we separate the diagonal and the non-diagonal 
contributions in Eq. llPQl) 

- 2 

Str {Q n kA) = Str2 (QnkA) + ^ Str (Q n kA) Str (Q p kA) (21) 

. n J n n^p 

For the diagonal contribution the integration over any Q m with m ^ n gives unity due to 
supersymmetry. Thus we obtain N equal terms which we denote by . For the same 
reason, the integration over Q m with m ^ n or m ^ p is trivial for the non-diagonal 
contribution. In this case we have iV 2 — N equal terms I^J and hence R^\u) can be 
written as 

2 



R<®( U ) =( — | Re 



r(0) 



T (o) 

L nd 



Nlf + (N 2 - iV)/ ( 
J dQ Str 2 (QkA)e^ st * AQ , 
J dQ Str (Qk A) StrA ^ 2 . 



(o) 

nd i 



(22) 
(23) 
(24) 



One can notice that the integral 1^ is exactly the same as the one, which determines 
i?2(^) in the Wigner-Dyson RMT. It can be easily calculated using the parametrization 



of Q given in Appendix A 



f) = l 6 + ^e 4 ^sin4fi, 



71 UU 



(25) 



This result corresponds to the Wigner-Dyson expression for H.2{u)) pi], which is not 
surprising, as we assume that the statistical properties of a single grain are described 
by the standard RMT. 

The calculation of the non-diagonal contribution is even simpler, as according 
to Eq. flA~5|) the expression for Str (QkA) does not contain the full set of Grassmann 
variables and hence the integral in Eq. (|24|) is anomalous and equals to a constant: 



I {0) - 16 

1 nd — i0 - 



(26) 
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Eq.f l22p shows that for Q ^ it is precisely this contribution that is dominant in 
the thermodynamic limit iV — > oo. Additionally we should take into account the self- 
correlation contribution contained in Eq.( l25l) . as it survives in the thermodynamic limit, 
provided that uN is kept constant. Thus we find that 

= 1 + -fiS(w), asA^oo. (27) 

There is a simple way to understand the above result without doing any calculations. 
Indeed, although the level correlations in each grain are non-trivial and described by 
the Wigner-Dyson formula, the levels in different grains are assumed to be completely 
uncorrelated in the zeroth order of the virial expansion. For this reason, when the 
number of grains A^ — > oo, all possible correlations except for the self-correlation vanish 
and (u(E - tu/2) v(E + u/2)) = (v(E - u/2))(u(E + u/2)), which leads to the trivial 
result (127]). provided that the density of states in an energy interval [E — u/2, E + u/2] 
can be approximated by a constant. 



4-2. Two-level correlation function in the first order approximation 

The first step in calculation of R^\uj) is to use Eq. (|2T|) to separate the diagonal and the 
non-diagonal contributions in Eq. (I20|) . Then one can realize that the diagonal terms give 
a non-zero result only, if n = m or n = k. All other terms are anomalous and vanish, 
as F m k = at Q m = Q^ = A. In a similar way, for the non-diagonal contribution we 
should keep only the terms with n = m and p = k or n = k and p = m. Thus we need 
to deal with integrals of two types: 

2 



4%') = i^i ^ 



i 



m>k 

J« = J dQdQ' Str 2 (QkA) (e~ bStTQQ ' - l) StoA (Q+Q') ; (29) 
J« = J dQdQ' Str (QkA) Str (Q'kA) (V^W _ ^ e ^strA(Q+Q') ; (30 ) 

where the dependence of 1^ and 1^] on m and k is hidden in the parameter b = 2U mk jt. 
Calculation of the above integrals is more challenging problem than the one we were 
faced in the previous section, as now we are not allowed to consider the limit u — > 0. 
The main steps of the calculations are outlined below and the details are presented in 



Appendix B 



We start with integration over the Grassmann variables and the phases in a similar 
way as it was done in the previous section. Then we are left with a sum of four-fold 
integrals over bosonic and fermionic variables Xg, X' B and Xp, X' F . The only term 
which couples the bosonic with the fermionic variables originates from the measure 
of integration. It can be represented as an integral over an auxiliary variable of the 
exponential function. In this way the bosonic and the fermionic integrals can be 
completely decoupled and considered separately. 
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The next step is to recall that we are interested in the strong coupling limit t — > oo 
or b 0. One can show that in this limit R^\oj) is given by the asymptotic series 



L 2 

fe>0 



■-2 

^m = EM b ^(^ ( 31 ) 



The higher order terms of the virial expansion R% (oo) have similar representations, but 
with k > i — 1. Thus only the first term of the series (l3Tj) is relevant in the first order 
approximation. For this reason, we are allowed to take the limit b — > keeping u/\/b 
constant. 

It turns out that the fermionic integrals can be easily calculated in this limit by 
expanding the corresponding integrands into the power series. All of them give just 
trivial constant contributions. The bosonic integrals are more tricky due to the fact 
that the bosonic variables are non-compact. Rescaling the integration variables and 
keeping only the leading terms of the expansion fl3T|) we arrive at two-fold integrals 
which can be evaluated exactly and gives the following expression for R^\uj): 



N 2 

m>k 



e 8U mk erfi 



(32) 



Adding this contribution to the formula for R^\uj) we obtain the final result for R 2 (u> 



*(«,) = -m + w E ^= e - erfi ^ j . (33) 

This is the main result of this section. Similarly to Eq. fTl5|) it can be used to 
compute i?2(w) for the a- model with an arbitrary coupling matrix U. Remarkably, 
the obtained expression coincides exactly with the corresponding result derived for the 
almost diagonal random matrices [8], if one identifies the coupling matrix U m k/t with a 
variance matrix n 2 u 2 (\H mk \ 2 \) /2(\H kk \ 2 ) . 



5. Critical power-law cr-model 

In this section we apply the general expressions for I q and Rziui) derived in the previous 
sections to the cr-model showing critical behavior. One of the most studied example 
of this kind is the model characterized by the power-law decay of the coupling matrix 
elements U mk oc \m — k\~ a . Such a model was derived in Ref.[H] as an effective field 
theory corresponding to the ensemble of the power-law banded random matrices. It 
is believed that it exhibits critical properties at all values of the coupling constant 
provided that a = 2. The mapping between the ensemble of the random matrices and 
the cx-model can be justified only in the large bandwidth (small coupling) limit. For 
this reason the critical cr-model was studied in detail in this limit [TU [151 02] ■ At the 
same time the model is well defined at arbitrary values of the coupling constant, in 
particular in the opposite regime of the strong coupling, where it can be explored by 
the method developed in this work. On the other hand, many analytical results are 
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available for the original random matrix model in the case of the small bandwidth limit 
p~5j \T7\ [HI [191 |20j [21] . This enables us to compare the predictions obtained for the 
microscopic model and for the cr-model in the situation, when it seems to be unfeasible 
to establish a direct relation between them. 

Specifically, we consider the cr-model (CQ) with the coupling matrix 



7T 



2 



-U mk = B 2 "-, r-. (34) 

The periodic form of the coupling matrix ( I34j) is more convenient for analytical 
calculations. 

In order to compute I q we need to plug the coupling matrix (|34|) into the general 
expression ( II 5p . Calculating the sum over m and n in the leading order in N we obtain 
the following expression 

h = r(g + 1) (l - B T ^~y^ In N + ■ ■ -V q>\. (35) 

Comparing it with the equation I q oc N^ dq ^ l " > defining the fractal dimensions d q we 
find the formula for d q in the lowest order in B: 

«. = B'-^p, (36) 

which is the same as the result that one would obtain for the power-law random model 
with (\H mm \ 2 ) = 1 and 

(\H mk \ 2 )= (^=) 2 f m^k. (37) 

We conclude that identifying the coupling matrix U m k/t in the cr-model with 
the expression 7r(|i^ m fc| 2 )/4 in the random matrix model, we obtain a one-to-one 
correspondence between two results. 

Next we calculate R2(w) for the same model. To this end we need to evaluate the 
double sum in Eq.( |32l) for U m k given by Eq.(j34j) in the limit iV — > oo. In this limit one 
can first introduce rescaled variables s = Nuu, x = sin(7r|m — k\/N), <fi = nm/N and 
then replace the sums by the corresponding integrals: 

(i), , 2s f 1 , , T in0 . x „ fsx\ 



Wis) = -l + —f dcP dx e-^erfi — . (38) 

Bit? Jo Jo y/l-x 2 V b J 

Integrating by parts over we are left with an integral over x, which can be evaluated 

exactly: 

2 

B£\s) = -e~^. (39) 
Finally, taking into account the contribution of R% (s) we obtain 

R 2 (s) = 8(s) + l-e-£. (40) 

Again we notice that this expression is in one-to-one correspondence with the random 
matrix result, provided that we use the same identification rule U m k/t — > n{\H mk \ 2 ) /4 
as before. 
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6. Conclusions 

The explicit calculations performed in the previous sections show that the first order 
results both for I q and i^t^O are equivalent to the corresponding expressions derived 
for random matrix models using the virial expansion method. The equivalence is 
very general and holds true for an arbitrary coupling matrix \J m k- In particular, we 
demonstrated that there is one-to-one correspondence between these two quantities 
calculated in the ensemble of the power-law banded random matrices at b <C 1 and their 
counterparts calculated in the cr-model. Such a correspondence is rather unexpected 
taking into account the fact, that the standard mapping of the the power-law banded 
random matrix ensemble onto the cr-model can be justified only in the opposite limit 
6» 1. 

A natural question, which appears in this context, is whether the observed 
equivalence can be extended to higher order terms and if it is the case, then what 
is the reason behind it? Below we show that under rather general assumptions both the 
cr-models and the random matrix models can be mapped onto the same effective field 
theory, so that such an equivalence is preserved at any order of the virial expansion. 

In Ref. [9 J it was noticed, that one can formulate a supersymmetric field-theoretical 
description of disordered quantum systems or random matrix models, which is similar 
to the non-linear cr-model, but it is exact and free from the saddle-point approximation. 
The form of the action of the new theory is similar to the action of the cr-model (JTJ, but 
the target space of the new supermatrix Q is different from the target space of Q. This 
approach was subsequently applied to the almost diagonal random matrices in Ref. [8]. 
In particular, it was shown in [5] that in the strong coupling limit the target space of Q 
can be reduced by the constraint StrQ = 0, provided that the large scale approximation 
can be used. The latter assumes that an exact form of the coupling matrix U}. m at small 
distances \k — m\ plays no role, and this condition is automatically satisfied for a large 
class of the coupling matrices in the thermodynamic limit iV — > oo. 

The constraint Str Q = implies that only two real and four Grassmann variables 
are required for parametrization of Q. Its block structure in the retarded-advanced 
representation is similar to the block structure of Q given by Eq. llA.lj) : 



U \ ( Drr d ra \ ( u- 1 



e H o v){dZ dZ){ o I. <«> 

where U and V are exactly the same matrices as those that involved in the Efetov's 
parametrization of Q and defined by Eq. (1A.2l) [9]. Each of the matrix D a p has only one 
non-zero matrix element: 

D RR = diag(A, 0), Dra = diag(-Ae"^, 0), 

Dar = diag(Ae^, 0), Drr = diag(-A, 0), 

Ae[0,oo), 0g[O,2tt). (42) 

Now let us return to the virial expansion of the cr-model. To be specific we focus on 
i?2(w), but our consideration can be applied for any other correlation function. The nth 
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term of the virial expansion involves integration over n supermatrices Qk and it can be 
written as the following asymptotic series: 

4%) = £(^)~ <t+ "~V(^<)- (43) 

fc>0 

Using the same arguments as those, which were employed in the discussion of the validity 
of the approximation Str Q = (see Ref.[8] for details), one can show that actually all 
the terms with k > 1 can be neglected, provided that the large scale approximation is 



valid. In |Appendix C| we discuss the validity of this approximation for the power-law 
cr-model. 

In order to find the first term of the series, we scale the variables A = Xs/Vt and 
Co = y/tco and take the limit t — > oo. In this limit the integrand preserves its form, but 
the matrices Qk should be replaced by 

1 

Qk = £P ~T + Qk I A is constant ■ ( 44 ) 

Looking at the parametrization of Q ( 1A.lHA.3l) we can immediately notice that Q can 



be obtained from Q by setting the variables A^ and pp equal to zero and by replacing 
Pb by \b- Identifying <p B with from Eq.( |42|) we conclude that Q = Q. The integrands 
for both models are the same and the measure can be written as 

dQ = dQ = 2j;7^ d P dp*dia dia*. (45) 

Since the above argument holds true at any order of the virial expansion we arrive at 
the conclusion that in the strong coupling limit and under the assumption of the large 
scale approximation both models can be mapped onto the same supersymmetric field 
theory. The action of this theory is the same as the one of the a-model and the field 
variables Qk are parametrized according to Eqs. fl41|42l) . In particular, it implies that 
all the results derived for the almost diagonal random matrix models are applicable for 
the corresponding cr-models in the strong coupling limit and vice versa. 
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Appendix A. Parametrization of supermatrix Q 

We employ the standard parametrization of Q introduced by Efetov [3] and adopt the 
notation used in |12j . In the unitary symmetry class (/3 = 2) Q has the following block 
structure in the retarded-advanced representation: 
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The matrices U and V contain only Grassmann variables p, p* and a, a*: 

uJ 1+ 0?' 2 , /, vJ 1 -^ 2 , * |.(A.2) 

I p* 1 + p p/2 I \ 10* 1 - a*a/2 ' 

The matrices D a/3 contain other four real variables X B , X F and 4> B , (fi F : 

d rr = -D AA = ( A Q B ^ V X B e [1, oo), X B e [-1, l] 

I -A*Be-* B \ _ / /i B e^ 

D/M " 1 ^e"*' j ' Afi ~ I 



B , F G [0,2tt), // b = yJX 2 B -l, p F = sJl~\ 2 F . (A.3) 
The measure in this parametrization is given by 

dQ = ^ B ^ F ^t, F dp dp*dia dia*. (A.4) 

Using the expressions for the matrices k = diag(l, — 1, 1, —1) and A = diag(l, 1, — 1, —1) 
and the above representation of Q, one can easily calculate Str(AQ) and Str(kAQ): 

Str(AQ) = 2(X B -X F ), 

Str(AAQ) = 2 [X B + X F + (X B - X F ) {pp* - aa*)} . (A.5) 

Appendix B. Two-level correlation function in the first order approximation 

According to Eq.f l29|) the diagonal contribution in the first order of the virial expansion 
can be written as 

/«=r(6)-T(0), (B.l) 
T{b) = J dQ Str 2 (QkA) e 2instvAQ Y(Q), fi = (B.2) 

Y(Q) = J d Q> e -bStrQQ>+2inStrAQ' (33) 

Changing the integration variable Q' by U~ 1 Q'U in the integral defining the function 
Y(Q), one can show that Y(Q) = Y(X B , X F ) [10j[T2]. Using this fact and the expressions 



for StrAQ and StrQkA from | Appendix A\ we obtain 

T = 4 r dX B f dX F - \ r= t^E^L f dp*dpdta*dtae 4in ^- x ^x 
Ji J-i (As - X F ) 2 J (27T) 2 J 

[{Xb + X F ) 2 + 2(A| - X 2 F )(pp* - 00*) - 2{X B - X F ) 2 pp*aa*] Y(X B , X F ). (B.4) 

The three terms in the brackets in the above expression generate three contributions to 
T. The first one T\ corresponding to [X B + A^) 2 is anomalous and therefore is given by 

ri = 16K(l,l). (B.5) 
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The second one corresponding to 2(A^ — X F )(pp* — aa*) is zero, as the bosonic part of 
the integral converges and the fermionic one is zero. Finally, the third one T 2 generated 
by the last term in the brackets has the form 

/oo /*1 
al\ B J d\ F e^ XB - XF) Y{\ B) \ F ). (B.6) 

Now we need to use the explicit form of the function Y, which can be found by 
performing integration in Eq.( lB.3j) . The result is available in Refs.|10[ IT2]: 

Y(X B ,X F ) = e - 2b ^ B -x F ) + 

d\' B f d\' F ^^£ e ^(A^A' F )- 2 6(A s A' fl -A,A^) L(A / B) ^ ^ ^ (R7) 
J 1 */ — 1 B F 

L(X' B , X' F , X B , X f ) = 2b 2 ((A' B A B + X' F X F )Io{2bp,' B p, B )I (2bfi' F p, F ))- 

[HBH'Bh(2bfj.' B HB)Io(2bn' F fj. F )) - p F p,' F h{2bp,' F p F )I {2bn' B p B ))}) , (B.8) 

where Jo and I\ are the Bessel functions. One can see that the above integral is zero at 
Xb = X F = 1 and hence Y(l, 1) = 1. Thus we find the answer for T\. 

T x = 16. (B.9) 

Another observation which can be made here is that the exponential function and the 
integral in Eq. (IB.7j) generate two different contributions in Eq. (IB.6j) . which we denote 
by T 2a and T 2 b respectively. The first one can be easily calculated: 

/oo /* 1 

dX B J dX F e ^{\ B -\ F ) e -2b{\ B -\ F ) = b ^ 

(2n + uy e ■ (R10) 

This exact result can be represented as a sum of two different contributions: the b 
independent part and the rest, which can be expanded into the power series: 

/(n,6) = /(n,o) + J2 (v^)' ' fk (4) • (b.ii) 

fc>o \ybj 

The explicit form of the b independent part is not important as it will be canceled out 
later according to Eq. ( TB.lj) . Moreover we are interested only in first term of the power 
series, as the higher order terms exceed the accuracy of the first order approximation of 
the virial expansion: 

isin(4Q) AiU b 
T 2a -— +2^ + .... (B.12) 

Calculation of T 2 b requires much more efforts, as T 2 b is given by a four-fold integral. 
However, one can notice that all the terms in Eq. (lB.7j) except from the first one are 
factorized into the fermionic and bosonic contributions. To exploit this observation we 
introduce the following function 

/OO /»1 /"OO /*1 

dX B / dX F / dX' B / dX' F e ^(\ B -\ F )+un'(\' B -x' F ) x 

e -2b(X B X> B -X F X> F ) L ^ x , Afl) Af); ( g 13) 



Virial expansion of the non-linear sigma model in the strong coupling limit 



15 



which can be used to find T 2 f, 

i r° , d „ ( % 



T » = uh dt an F { Q - n ' + t 



(B.14) 



It has an advantage that now the fermionic and the bosonic integrals do not mix up: 

4 />oo pi poo pi 

d y e 4in(x B -x F )+4iQ'(\' B -\' F ) x 



/oo pi poo p 

dX B J d\ F J dX'sj 



-1 

e -^ B x B -^' F ) fa{XBj \' B ) 9a (X F , X' F ), (B.15) 
fi = X' B X B I (2bn' B n B ), g 1 = I (2bfj,' F fj, F ), 
/a = I (2bu B u B ), 92 = X' F X F I (2bu F u F ), 
/s = -fi' B fi B Io(^bfi' B /j, B ), g 3 = I (2b[i' F fi F ), 

fi = I Q (2byl B n B ), #4 = n' F n F h{2byl F n F ). (B.16) 

Thus the problem is reduced now to the calculation of the two-fold bosonic and fermionic 
integrals, which we denote by A a and B a respectively: 

4 

F(Q, Q') = 16b 2 KB a . (B.17) 

a=l 

Again we are interested only in the leading term of the expansion (IB.llj) and therefore 
we are allowed to consider the limit b — » keeping the ratios Q/y/b and 0! / Vb constant. 
In this limit the fermionic integrals can be easily evaluated by expanding the integrands 
into the power series and the results read 

B 1 = B 3 =4 + 0(b), B 2 = B 4 = 0{b). (B.18) 

The bosonic integrals are less trivial due to the non-compact character of the bosonic 
variables. In order to expand the integrals into the power series we introduce new 
variables u = 2bX B X' B , v = X B /X' B , r\ = Q/vb and rj' = Q'/y/b. Then keeping only the 
leading order in b term in the integrands we obtain 



-1 POO POO 

Ax = 4= du cfo VW¥+Vv^)e-%( M ). 
8b 2 Jo Jo v 

A 3 =-4- du dv-e^VW+^)e- u h(u). (B.19) 
8b 2 Jo Jo v 

In a similar way, one can show that A^ converges while A 2 logarithmically diverges as 
b — > 0. According to Eq. (IB.18l) this implies that their contributions to F(Q,Q') can 
be neglected. Substituting the results of Eqs. flBT9l IB~T8l IB~T7j) into Eq. flBTiD and 
integrating over t and v we find 

T 2b = -32 J du uK 2 (^-SiVL^^j e- u \I x {u) - I (u)] , (B.20) 

where K 2 is the modified Bessel function. Collecting the results of Eqs. (IB.9llB.12IIB.14l) 
we obtain 

zsin(4f2) A - a b 
T(b) = 16 + —^r 1 ^ + 2- - 

32 jf du uK 2 ^~ 8i ^\J^j e ~ U t J i( M ) - J o(«)] • (B.21) 
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In order to compute the non-diagonal contribution Eq. (l30p exactly the same 
approach can be used. Skipping the details of the calculations, we present here the 
final result for S(b), which plays the same role as T(b) for the diagonal contribution: 



S(b) = 16-32 y du uK I -8iQ^ )? 

It turns out that both integrals in Eq. (1B.21I) and Eq.f lB.22j) can be computed analytically 
and the real part of their sum is given by 




[h{u) -/„(«)] 



(B.22) 



/■oo 

Re / duu 
Jo 



K 
b 

160 2 




Sii'UI^ )+K 2 



-8iO 



1 JnO _4o£ „ (20 
=- e b erh —= 

2 Vb \y/b 




e- u [h(u) - I (u)} 



(B.23) 



where erfi(,2) = (2/^/tt) f*e x2 dx. Therefore the final expression for the sum of the 
diagonal and non-diagonal contributions reads 



nd i 



-16 



1 - 



y/nO _4n£ n (20 
— — e b erh —= 



(B.24) 



Substituting this expression into Eq.( l28|) and recalling the definition of b and O we arrive 
at Eq.(|32|). 



Appendix C. Validity of the large scale approximation for the power-law 
a-model 



In this Appendix we show, on the example of the power-law cr-model, why the term of 
the order of 1/t in the expansion of 1^ can be neglected in comparison with the term 
of the same order originating from the expansion of Ig . 

It follows from Eqs. ( fl0lfT3|) that the next term in Eq. ( fl5|) has the following structure: 



u„ 



t 



(C.l) 



where G(q) is some function of q, which does not depend on N. Substituting the 
expression for U mn fJMjl into this formula we obtain that in the thermodynamic limit 
iV — > oo the result is CG(q)B 2 with some constant C. 

On the other hand the leading order in B term for I q will be the same as the one 
calculated for the almost diagonal random matrices and it is given according to Ref . [18] 
by the following expression 

B 2 {G x {q) In 2 N + G 2 {q) In N + . . .), (C.2) 

with some N- independent functions Gi(q) and G 2 (q)- One can see that although 
this term has the same order in B it will be dominant in the thermodynamic limit 
iV — > oo. The reason for a different iV-dependence appearing here is related to the 
fact, that in Iq there is an additional summation over the indices of the coupling 



matrix Emnpy t \/"f ■ ^ ne importance of such an additional summation for the 
justification of the large scale approximation was emphasized in Ref. [5] 
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